Creation of a quantum oscillator by classical control 
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As a pure quantum state is being approached via linear feedback, and the occupation number ap- 
proaches and eventually goes below unity, optimal control becomes crucial. We obtain theoretically 
the optimal feedback controller that minimizes the uncertainty for a general linear measurement 
process, and show that even in the absence of classical noise, a pure quantum state is not always 
achievable via feedback. For Markovian measurements, the deviation from minimum Heisenberg 
Uncertainty is found to be closely related to the extent to which the device beats the free-mass 
Standard Quantum Limit for force measurement. We then specialize to optical Markovian mea- 
surements, and demonstrate that a slight modification to the usual input-output scheme — either 
injecting frequency independent squeezed vacuum or making a homodyne detection at a non-phase 
quadrature — allows controlled states of kilogram-scale mirrors in future LIGO interferometers to 
reach occupation numbers significantly below unity. 



Motivation. The detectors of the Laser Interferome- 
ter Gravitational-wave Observatory (LIGO) [1] arc cur- 
rently operating at a factor of 10 above the Standard 
Quantum Limit (SQL) [2] at closest approach, limited 
by classical noise below about 100 Hz, and by quan- 
tum shot noise above 100 Hz [3]. This low noise perfor- 
mance allows for probing and manipulating LIGO mir- 
rors (10 kg each and suspended as pendulums with reso- 
nant frequency ujp = 27r x 0.7 Hz and typical damping rate 
7p ~ 10~® sec^^) close to scales set by the Heisenberg Un- 
certainty Principle. Recently, electronic feedback control 
was used to shift the resonant frequency of the pendulum 
mode up to ^ 140 Hz and damp it, leading to an effective 
occupation number of 234 for the kg-scale oscillator [4]. 
Semiclassical calculations estimate an effective occupa- 
tion number of TVoff « QcsSxi^cs)/SsQL{^cs), where 
Sx is the detector's position-referred noise spectral den- 
sity, S'sQL = 2h/{mVL^) is the free-mass SQL for position 
measurement, Vl^s is the eigenfrequency of the controlled 
oscillator, and Qcff is its quality factor. Extrapolation to 
A'eff < 1 favors Qofi <C 1; yet this is exactly where the 
approximation fails, and a more careful treatment of op- 
timal control is required. Moreover, we must answer the 
following question: does the strong continuous position 
measurement required for making the stiff electro- optical 
potential always produce significant decoherence of the os- 
cillator's quantum state? 

The cold damping technique used in Ref. [4] was first 
proposed by Mancini et al. [5] and demonstrated ex- 
perimentally by Cohadon et al. [6]. Subsequent exper- 
iments have used cold damping to cool mechanical os- 
cillators with the goal of reaching the quantum ground 
state [4, 6, 7, 8, 9, 10, 11, 12], typically using simple 
(proportional plus derivative) filters. Little attention has 



been paid to the question of whether these filters are 
optimal. Furthermore, feedback control was used only 
for damping in these experiments, and not for shifting 
the resonant frequency of the oscillator, as in the LIGO 
experiment [4]. In this Letter, we obtain the optimal 
state-preparation control scheme for a general linear mea- 
surement process (possibly non-Markovian) , and study 
prospects of quantum-state preparation for kg-scale test 
masses in future LIGO detectors via feedback control. 
Our general theory also applies to other mechanical struc- 
tures [6, 7, 8, 9, 10, 11, 12]. 

It turns out that occupation number is not necessarily 
a good measure of "quantum-ncss" : squeezed states, for 
example, can have high occupation number, yet they are 
more "quantum" than a vacuum state. Moreover, the 
definition of occupation number requires a well-defined 
real-valued eigenfrequency, which can be ambiguous for 
two reasons: (i) the controller can modify the oscillator's 
original eigenfrequency ujp] and (ii) for oscillators with a 
finite quality factor Qoff, the choice for an effective real 
eigenfrequency Q.cs would be ambiguous by ~ f^cff/Qcff- 
Instead we use the purity, defined as 

U = ^VxxVpp - > fi/2. (1) 

where Vxx, Vpp and Vxp are uncertainties of oscilla- 
tor's position and momentum along with their cross- 
correlation. For steady states, which have Vxp = 0, U 
can be converted to an effective occupation number, 

iVeff = U/h - 1/2 , (2) 

which is the minimum occupation number one could ob- 
tain when the same mirror state is put into a quadratic 
potential well with an arbitrary eigenfrequency Q: 

N,s + 1/2 = min [{Vpp/{2m) + mil^Vxx/2) /(hn)] , (3) 
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achieved at = \/Vpp/ (m'^Vxx)- Since ft may be 
very different from the oscillator eigenfrequency, the re- 
sulting quantum state tends to be position squeezed if 
51 > Wp, and momentum squeezed ii il < u>p [17]. iYeg 
also determines the von Neumann entropy of the state: 

S = (iVeff + 1) l0g(A^eff + 1) - iVeff logTVeff. 

General Optimal Controller. A block diagram of 
the entire measurement-control system is shown in Fig. 1, 
where x is the position of the oscillator, y is the output 
field we measure, H is the measurement transfer function, 
Rxx = — 1/[(^^ — ujp + ijp){fl + ujp + i^p)] is the response 
function of the oscillator (with 7p = uJp/Qp its relaxation 
rate; the tilde denotes a frequency-domain quantity), F 
and Z are force and sensing noises, G is a possible clas- 
sical force acting on the oscillator, and C is the feedback 
kernel. We write the closed-loop position and momentum 
of the oscillator as 

ictrl = So - ^ctriyo , Petri = PO + iflKctilVO (4) 

where xq, po are the open- loop evolution of the oscillator 
position and momentum (we set its mass equal to 1, and 
use p = —iflx), yo is the open- loop out-going field, and 
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(5) 



The closed-loop dynamics is stable and the feed- 
back C proper, if and only if (i) A'ctri is causal (i.e., 
no poles in the upper-half complex plane); and (ii) 
limo-,00 ftKcti\{ft) = 0. The closed-loop response of the 
oscillator's position to external force is given by 
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In Eq. (4) , closed- loop quantities are viewed as subtract- 
ing the open-loop readout field from open- loop quantities. 
On the other hand, causal Wiener filters Ka (a = x,p) 
can be constructed based on the cross spectral density be- 
tween xq, Po a-nd yo, to yield the best (least-mcan-square) 
estimates of xq and po based on past measurement of yo 
[16]. In terms of these filters, we have 
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where Syy is the single-sided spectral density of yo [18]. 
Minimizing over all stable closed-loop systems, we obtain 



[/opt ^ mm y^^^H^ = ^VsJ^p + V^p , (8) 
which is achieved by a unique controller with 
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Here 4>+^{n) is a causal whitening filter for the out- 
put field yo, with = Syy and both 0+ and 01]^^ 
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FIG. 1: Block diagram of the feedback control system. 

analytic in the upper-half complex plane, while Gx = 
[Sxy/4>+] + y where Sxy is the cross spectral density be- 
tween xo and yo, and [i^]+ denotes extracting the causal 
part of F while F* stands for the complex conjugate 
of F. More specifically, under the scaling transform 
{x',p'} = {x/ ^/V^,p/ ^yV^}, the error ellipse of the 
optimally controlled state becomes a circle with V^}"^] = 
Vp,^p} = 1-1- V^,p,, which is in turn equal to the larger 
eigenvalue of the re-scaled conditional covariance matrix. 
Thus the error ellipse of the optimally controlled state is 
the one with the minimum area among those that (i) to- 
tally encompass the conditional-state error ellipse; and 
(ii) are consistent with V^p '^ — 0. This means the con- 
trolled state is always a mixed state, unless V^p ~ [19]. 

General Markovian measurements. We consider 
the open-loop system (Fig. 1) 
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— Rxx{F + G) 



(10) 



where we have set H ^ I (without loss of generality), 
and Z and F are now characterized by real and constant 
(single-sided) cross spectral densities Szz, Szf and Sff 
satisfying the Heisenberg Uncertainty Relation [2] 
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where /i > 1 measures the purity of the measurement pro- 
cess (with /i = 1 corresponding to quantum-limited mea- 
surement). This describes measurement systems with 
white sensing and force noises, e.g., measurement of mir- 
ror location using a broadband Fabry-Perot cavity, with 
frequency independent input squeezing and homodyne 
detection. We recast Syy into a causally factorized form: 
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= SzzLL*/{PP*). 

Here P = -1/Rxx, LL* = 0.'^ - 2AujIQ? 

1 Szf „2 1,2 Szf 



Sff\Rx 



A=l 



i^p Szz 



1 



t^p Szz 



1 Sff 
^p Szz 



(12) 



(13) 



L is defined here in such a way that it only has roots in 
the upper-half complex plane [20]. Using the definition 
for Ka and Eqs. (12), (13), the conditional covariance 
matrix can be obtained: 
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The conditional purity is given by 



(15) 



which is identical to the "purity" of the measurement 
process; see Eq. (11). In the absence of classical noise, 
the conditional quantum state of the oscillator is always 
pure. With Eq. (8), we obtain 



U,u-i/{fi/2) = ^i{^l-A/B + V2)/Vl + A/B (16) 

which is achieved by the unique optimal filter [cf. (5) and 
(9)] with associated closed-loop response function 

c^Co{n-Ci)/{n-C2), (i7) 
Rf, = -{n - n4)/[{n - - n2)in - n^)], (is) 

where Co = -(c^^ + ^4^13), Ci = {n^ + ojIVLa)I{uoI + 
U.^^^), C2 ~ and f2i,2 are roots of Q, namely 



f}i,2 = ±uJp^/{B + A) 12 - iUp^{B - A)/2 (19) 
while 573^4 are purely imaginary: 



-iVBujp , ili = -i[VB + ^2{B - A)]ujp . (20) 



We note that: (i) the poles of the closed-loop dynamics, 
ill 2, a-re identical to the zeros of Syy, i.e., the optimal 
controller "finds" the frequency of maximal sensitivity, 
and shifts the oscillator's eigenfrcquency there; (ii) the 
51 that achieves A'off in Eq. (3) and motivates N^a as 
an occupation number in a harmonic potential, is equal 
to |Sli,2|! the modulus of the closed-loop poles; (iii) the 
optimal controller (17) is in fact proportional feedback 
plus constant damping and simple band limiting (which is 
required for V^'^p'' to be finite), justifying previous choices 
[4, 6, 7, 8, 9, 10, 11, 12]; and (iv) a pure state is only 
strictly achievable when ^ « 1 (i.e., in the absence of 
classical noise) and Ak, B {va general |^| < B). 

Having A k B corresponds to a high quality factor for 
the closed-loop dynamics [Cf. Eq. (19)], with 



= ^bTa/{2Vb^). 



(21) 



This is consistent with our understanding that a low-Q 
oscillator cannot have a pure quantum state due to the 
Fluctuation Dissipation Theorem [2]. Moreover, A B 
also corresponds to Syy having a very small minimum [in 
the limit oi A = B, Syy reaches 0; cf. Eq. (12)]. Let us 
consider the force (G)-refcrred noise spectrum, 



5, 



G 



^yy 1 1 I 



— SzzLL* ^ 



(22) 



and compare it with the free-mass SQL for force detection 
Sq^^ = 2n^h. Taking the minimum over all frequencies, 
we obtain the factor by which the SQL is beaten: 
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FIG. 2: Left panel: N^s as a function of measurement 
strength Qg/tOp. Curves (a-e) represent A'efi at different tem- 
peratures 6: a) 0.1, b) 0.5, c) 1 (critical), d) 2, e) 10. Right 
panel: Minimal A^efi as a function of factor 77^1 for 10 dB 
squeezing (blue) and vacuum input (black). 



which leads to [cf. Eq. (16)] 
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>V^ + l, iVeff > 77V2. (24) 



This means an oscillator under measurement can only 
be converted into a quantum oscillator via control if it 
can beat the free-mass SQL in a force measurement, in 
which case the optimally-controlled closed-loop quality 
factor Qcff would also far exceed unity. 

Position measurement with light. In the realis- 
tic case with suspension and internal thermal noises and 
optical loss, we have (as in Ref. [14]) 

X = Rr,^,[aai + (,F + G] (25a) 
y = aiSUKf) + cos(l)[a2 + a/h{x + Cx)] + V^n (25b) 



or in the notation of Eq. (10): 



F = aai +^_F, Z = 



ai sin (j> + a2 cos (p + y/en 
a/hcoscf) 



,(26) 



where ai.2 are the input quadrature fields, is the read- 
out phase (0 for phase quadrature and tt/2 for ampli- 
tude quadrature), a is the measurement strength (a = 
i^^hujolc/ (tc^) for a Michelson interferometer with arm 
cavities, with wq the carrier frequency, Ic the circulating 
power in the arms, r the input-mirror power transmis- 
sivity, c the speed of light), e is the optical loss, and n 
is the vacuum noise, and are the classical sensing 
and force noises, respectively, whose spectra cross the 
position- and force-SQL at frequencies Qq/Cx and ^qCf, 
respectively, where Jig = a/^/h is the characteristic mea- 
surement frequency (as defined in Ref. [14]). 

Viscosity noise alone and phase-quadrature 
readout. Here we consider the cold damped systems 
of Refs. [4, 6, 7, 8, 9, 10, 11, 12], i.e. e = = = 0, 
LUp ^ 0, Cf = 47p/cBT/(fl^ft). Since a free mass under 
such measurements never beats the free- mass SQL, ac- 
cording to Eq. (24), it will always have iVoff > 1/2. In 
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order to have a lower N^s, we have to rely on beating the 
free-mass SQL around the oscillator's original resonant 
frequency, which is only possible when kgT/ (hujpQp) < 
1/2. Analytic results reveal a phase-transition-like situ- 
ation: below a critical temperature, with 

e = T/T, < 1, T, = ha;pQp/{2V2kB) , (27) 
the minimum occupation number 



NoptiO) = 2-^/^[V2-6^ + ^29V2-e^ + 9-V2] (28) 
Nopt 2-3/46)1/2 as 6* 0, can be achieved with 



UJr, 



91/2(2-612)3/4 

V2~¥-9' 



-|^(^.)V4 



[0 ^ 0), 



while for T > Tc, a temperature- independent minimum 
of Ncs = is achieved with infinite measurement 

strength as indicated in Fig. 2 (left panel). 

For a LIGO detector, T > so iVoff > 1/V2 even if 
viscous damping alone is considered. For the systems of 
Refs. [6, 7, 8, 9], prospects for surpassing Tc — 17 K x 
(Qp/lO^) X [ujp/{2Tr X IMHz)] are much more promising. 

Prospects for LIGO {ojp = 0). In order to 
evade the limitation of Nes > l/%/2, we consider an 
arbitrary readout quadrature with squeezed light input. 
Given a classical noise budget, with vacuum input or 
fixed squeezing factor (10 dB), we optimize flq, 4> and 
the squeezing angle for a minimum iVefj. In the right 
panel of Fig. 2 , we plot A'cff as a function of the factor 
= [S'=V5'^^^]min = KfCx by which the total classi- 
cal noise beats the SQL [14] , while fixing the optical loss 
e = 0.01. The input squeezing appears to be not so cru- 
cial and almost the same results can be achieved without 
it. A detailed optimization will be presented elsewhere. 

Summary In this Letter, we developed the gen- 
eral theory for optimal state-preparation via linear feed- 
back control. The optimally controlled x-p error ellipse, 
achievable using the unique optimal controller, is the one 
with minimum area that still maintains Vxp = and en- 
compasses the conditional-state error ellipse. For a gen- 
eral Markovian measurement process, the conditional- 
state purity equals that of the measurement [cf. Eq. (15)], 
and the absence of classical noise guarantees a pure con- 
ditional state; yet a nearly pure controlled state requires 
additionally that correlations exist between sensing and 
force noises, in such a way that the device beats the free- 
mass SQL significantly in a force measurement. In this 
case, the optimal controller creates an oscillator with an 
eigenfrequency fieff around the most sensitive frequency 
of the device, with quality factor Qeff ^ 1 and an effec- 
tive occupation number A'eff <SC 1. This A^eff is also mean- 
ingfully measured against a harmonic potential with real- 
valued eigenfrequency w f^off- Furthermore, restrict- 
ing to conventional measurements applied to a viscously 



damped oscillator, we found a critical temperature T^, 
above which the oscillator is limited by A"eff > l/%/2, 
while below which, A"eff approaches as T/Tc 0. Fi- 
nally, for LIGO, the answer to the question raised at the 
beginning of this Letter is "yes" for conventional phase- 
quadrature readout, where the need to up-shift the os- 
cillator's eigenfrequency does limit N^s above l/-\/2, and 
"no" for an optimal choice of measurement strength and 
(non-phase) readout quadrature (possibly in combination 
with input squeezing), where the sub-SQL classical noise 
budget of future detectors will allow an optimally con- 
trolled kg-scale oscillator to achieve A"off ^ 1. 
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